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Abstract 

We study the heavy quark evolution in a quark-gluon plasma medium within the framework of 
Langevin equation coupled to a (2+l)-dimensional viscous hydrodynamic model. We modify 
the current Langevin algorithm such that apart from quasi-elastic scattering, medium-induced 
radiative energy loss is incorporated as well by treating gluon radiation as an extra force term. 
We find a significant effect of gluon radiation on heavy quark energy loss at LHC energies. 
Our calculation provides a good description of the D meson suppression measured by ALICE 
experiment, and makes a prediction for B meson suppression and flow. 



1. Introduction 

A highly excited deconfined state of QCD matter, the strongly interacting quark-gluon plasma 
(sQGP), is believed to be produced in relativistic heavy-ion collisions at RHIC and LHC. It dis- 
plays properties similar to a nearly perfect fluid and has been successfully described by hydro- 
dynamic models. Among various probes of the medium properties, heavy flavor has generated 
significant interest in the past decade, since heavy flavor high-p^ suppression and elliptic flow 
have been revealed as significant as for light quarks in spite of the large mass HHHI^l], indicating 
a stronger coupling strength between heavy quark and QGP than expected. Therefore, it appears 
crucial to explore how heavy quarks interact with such a hot and dense medium. 

Between the two energy loss mechanisms of heavy quark, gluon radiation is usually consid- 
ered negligible compared with quasi-elastic scattering, as long as the heavy quark energy is not 
very large, because of the "dead-cone effect" [4]. In the limit of multiple scatterings where the 
momentum transfer during each interaction is small, collisional energy loss inside a thermalized 
medium can be well described by the Langevin equation . In this framework, (60 provides 
a nice description of heavy flavor suppression and elliptic flow measured by RHIC. However, as 
we extend such study to the LHC energy level, even heavy quarks become ultra-relativistic, and 
therefore, it may no longer be reasonable to ignore the radiative energy loss due to the "dead-cone 
effect". 

In this paper, we shall modify the current Langevin approach such that gluon radiation is also 
incorporated as an extra force term, whose momentum space distribution will be governed by 
the Higher- Twist calculation 0, S Ht] . Within this improved approach, the competition between 
the two energy loss mechanisms will be clearly displayed as the energy increases. Combined 
with a (2+l)-dimensional hydrodynamic simulation of the QGP medium |T(i IT, 12], we find 
our calculation of D meson suppression to be consistent with the LHC data. 
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2. Methodology 



We modify the Langevin equation as follow: 



dp 
dt 



— = -t]d(p)P + £ + f s - (1) 



The first two terms on the right are the drag force and the thermal random force from the original 
Langevin equation for Brownian motion, and the third term f g = -dp g /dt is introduced for the 
force exerted on the heavy quark by gluon radiation. The probability of gluon radiation and 
the momentum distribution of the radiated gluon (Ap g ) during each time interval Af is sampled 
according to Eq.(0 taken from the Higher- Twist calculation 0. fi S : 

dN g _2a^_ p( a 2 lt zIi \( *j \\ 



dxdk\dt 7i k ± \ 2tj j\k\ + x 2 M 2 

in which q is the gluon transport coefficient, k ± is the transverse momentum of the radiated gluon, 
and x is the ratio between the gluon energy co and the heavy quark energy E. Additionally, t/ 
represents the gluon formation time and P(x) is the splitting function. Note that the last bracketed 
term in Eq.© is the "dead cone factor" for a heavy quark with mass M. 

With the assumption that the interaction during each scattering is small, the fluctuation- 
dissipation relation between the drag and the thermal force still holds - rjoip) = k/(2TE), where 
k is the momentum space diffusion coefficient defined in (^'(f)f (O) — K5'^5{t - t'). Meanwhile, 
we set the lower limit of the gluon energy to be o>q = 7tT, which is the balancing point between 
gluon emission and absorption. These assumptions guarantee the equilibrium of heavy quarks 
after a sufficiently long time of evolution in the medium. The different transport coefficients are 
related via D = 2T 2 /k and q = 2kCa/Cf, where D is the spatial diffusion coefficient of the heavy 
quark in the QGP. 

We shall use our modified Langevin equation to simulate the heavy quark evolution inside a 
QGP. The QGP medium is generated with a (2+ 1 )-dimensional hydrodynamic model lllOl 
with the MC-Glauber initialization. And the coupling of the Langevin approach to the QGP is 
in the same way as our previous study lfl3l [l4ll . Our heavy quarks are initialized with the MC- 
Glauber model for the position space and a leading-order pQCD calculation for the momentum 
space. After traversing the medium, they are fragmented to heavy mesons via Pythia 6.4 lfl5ll . 



3. Heavy flavor evolution, suppression and flow 

We start with the evolution of the charm quark energy spectrum due to the different energy 
loss mechanisms. Here, the charm quarks are all initialized with 15 GeV energy before traveling 
through an infinite medium with a fixed temperature of 300 MeV. The spatial diffusion coefficient 
in our study is set at D = 6/(2nT), corresponding to a gluon transport coefficient q around 
1.3 GeV 2 /fm at 300 MeV. As shown by FigQ] before 2 fm/c, collisional energy loss dominates 
the charm quark evolution. However, after 2 fm/c, gluon radiation starts to dominate. 

Figf2] displays the energy loss of charm quarks after propagating through a realistic QGP 
medium produced in Pb-Pb collisions at the LHC energy level. The x-axis represents the initial 
energy of charm quarks and y represents the energy loss. One observes that while collisional 
energy loss dominates the low energy region, gluon radiation dominates the high energy region. 
The crossing point is around 6 GeV, indicating that the collisional energy loss alone may describe 

2 



0.8 
0.6 

i M 

i 

0.2 



0.6 

1 0.4 

0.2 


0.6 

§ 0.4 

s 

0.2 



■ 1 fm/c 

. 2 fm/c 

. 3 fm/c 

. ■ — ■ 4 fm/c 


1 1 1 1 1 1 

Collisional / 

-Mr, 
/■< h 
/•'/ / • 

, . , y,jJ. 


\ 
\ 

A - 




Radiative 

( 

, ; * >s /' .1 




Total 

„ i r. ■ ■< - - -i ■ • ■ 


i 

.i 
\ 



3 6 9 12 15 

E (GeV) 



Figure 1 : (Color online) A comparison of the evolution 
of the charm quark energy distribution in a static medium 
between collisional, radiative and total energy loss. 
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Figure 2: (Color online) A comparison of the charm 
quark energy loss after they traverse the QGP medium 
created in central Pb-Pb collisions with V* = 2.76 TeV 
between collisional energy loss only, radiative only and 
the total effect. The crossing point between the colli- 
sional dominating region and the radiative dominating 
region is around 6 GeV. 



the heavy flavor observed by the current RHIC measurement well, but will become insufficient 
when we extend to LHC energies. 

Fig[3] shows our calculation of D meson suppression compared with the LHC data. Neither 
collisional nor radiative energy loss alone are able to describe the data. However, our combina- 
tion of the two mechanisms provide a good description of D meson Raa measured by ALICE. We 
also present our calculation of D meson \>2 in FiglU and our prediction for B meson suppression 
and flow in Fig[5]and Fig|6] B mesons display similar behaviors as D mesons except that the 
crossing point between collisional and radiative dominating regions becomes much higher due 
to the larger mass of the bottom quark vs. the charm quark. 
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Figure 3: (Color online) A comparison of D meson R^a 
between different mechanisms and the LHC data. 



Figure 4: (Color online) A comparison of D meson v% be- 
tween between different mechanisms and the LHC data. 
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Figure 5: (Color online) A prediction of B meson Raa 
with different energy loss mechanisms. 



Figure 6: (Color online) A prediction of B meson V2 with 
different energy loss mechanisms. 



4. Summary 

We have studied heavy quark evolution and energy loss inside a QGP medium in the frame- 
work of a Langevin equation. We have improved the current Langevin approach such that not 
only collisional, but also radiative energy loss is able to be incorporated by treating gluon radi- 
ation as an extra force term. Our calculation reveals a significant effect of the medium-induced 
gluon radiation on the heavy quark energy loss at LHC energies, which should no longer be ne- 
glected due to the "dead-cone effect". And our combination of the two energy loss mechanisms 
provides a good description of the D meson suppression measured by the ALICE Collaboration. 
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